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We study the impact of cosmological parameters' uncertainties on estimates of the primordial NG 
parameter /nl in local and equilateral models of non-Gaussianity. We show that propagating 
these errors increases the /nl la uncertainty by a term 5/n£''V/nl''' - 16% for WMAP and 
^/nl^'V/nl"' — 5% for Planck in the local case, whereas for equilateral configurations the cor- 
rection term are 5/nl"'V/nl"'' — ^^'^ "^/nl^V/nl"' — respectively. If we assume for /nl''' 
a central value (/nl^') — 60, according to recent WMAP 5-years estimates, we obtain for Planck a 
final correction A/^^^' ~ 3. Although not dramatic, this correction is at the level of the expected 
estimator uncertainty for Planck, and should then be taken into account when quoting the signifi- 
cance of an eventual future detection. In current estimates of /nl the cosmological parameters are 
held fixed at their best-fit values. We finally note that the impact of uncertainties in the cosmolog- 
ical parameters on the final /nl error bar would become totally negligible if the parameters were 
allowed to vary in the analysis and then marginalized over. 



I. INTRODUCTION 

Research of primordial non-Gaussianity in the Cosmic Microwave Background (CMB) is an important field of cos- 
mology today. In a recent work, Yadav and Wandelt [jy] claim a 3a detection of a large NG signal in the WMAP 
3-years data. The following WMAP 5-years analysis produced similar results but showed a reduction of statistical 
significance from about 3 to 2a [2]. Yadav and Wandelt estimate, using WMAP 3-years data, is 

27 < fj^l''^ < 147 (95%c.l.) . (1) 

The WMAP 5-years estimate is: 

-9</ilj°r'<lll (95%c.l.), (2) 

where /nl is a dimensionless parameter defining the strength of primordial non-Gaussianity (NG) [3i], and the super- 
script "local" indicates that we are considering a primordial NG curvature perturbation of the form 

<i>(x) = <i>L(x) + /^°r' (*!(x) - m^))) , (3) 

where $(x) is the total gravitational potential, <i>L(x) is the gravitational potential computed at the linear level and 
(• • •) stands for the ensemble average. The reason for calling this a local NG model relies in the fact that the NG 
part of the primordial curvature perturbation is a local functional of the Gaussian part. The local shape of NG arises 
from standard single-field slow-roll inflation [4] as well as from alternative inflationary scenarios for the generation of 
primordial perturbations, like the curvaton or inhomogenous (pre)reheating models @, or even from alternatives 
to inflation, such as ekpyrotic and cyclic models fi\. Other models, such as DBI inflation fsj and ghost inflation [oj, 
predict a different kind of primordial NG, called "equilateral", because the three point function for this kind of NG 
is peaked on equilateral configurations, in which the lengths of the three wavevectors forming a triangle in Fourier 
space are equal |T^. This second form of NG is characterised by the parameter f^^^, which defines the amplitude of 
the equilateral triangles. In the following we will sometimes write the amplitude of NG simply as /nl, without any 
superscript. When we do so, we mean that our conclusions apply to both the local and the equilateral case, with no 
need for distinction. 

Standard single field inflation predicts /nl''' ^ 10^^ at the end of inflation 0] (and therefore a flnal value /nl^' ^ 
unity after general relativistic second-order perturbation effetcs are taken into account [HI). It is thus clear that 
large central values of /nl^', like those obtained in the above mentioned analyses, are going to rule out the simplest 
scenarios of inflation as viable models of the Early Universe. On the other hand, the low statistical significance of 
the final WMAP 5-years result makes any conclusion premature at this stage. With its high angular resolution and 
sensitivity Planck will allow to significantly improve the statistical estimate of /nl'^', reducing the final error bars 
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from the present A/^^J^^^ ~ 30 to a final value of Afj^l""^ ~ 5 fl^ thus allowing a many a detection of non-Gaussianity 
if the present large central values of were to be confirmed. An eventual detection of a large /nl by Planck 

would not however automatically imply that the observed non-Gaussianity is primordial in origin. A number of effects 
can produce a spurious NG signal that can bias the final estimate. The most relevant examples of NG contaminants 
are probably given by diffuse foregrounds emission, unresolved point sources contamination, NG noise. Both the 
analyses by Yadav and Wandelt and by the WMAP team consider all these effects and conclude that they do not 
significantly affect the fj^f^^^ measurement from WMAP data. The picture gets however more complicated when we 
consider future Planck data. It has been recently shown that several effects that are not important in the analysis of 
WMAP data become no longer negligible with Planck. For example, Serra and Gooray [l^ studied NG contributions 
arising from several secondary sources, and concluded that effects such as the cross-correlations SZ-lensing and ISW- 
lensing produce a bias in the estimate of /nl*^' which is at the level of the expected estimator variance at Planck 
angular resolution. Analogous conclusions have been reached by Babich and Pierpaoli fT2| for the cross correlations 
of density and lensing magnification of radio and SZ point sources with the ISW effect. Note that all these effects 
are unimportant for present analyses both because they involve higher multipoles than those reached at the WMAP 
angular resolution, and because they produce a bias Afj^f^'^^ ^ 1, thus much smaller than the WMAP sensitivity 
to /nl'*': which is A/nl**' ^ 30. However, the much higher angular resolution achieved by Planck and an expected 
predicted sensitivity on /nl**' given by Af^l""^ - 5 (and A/^lff^' - 3 ff polarisation data are included in the analysis) 
will make the above mentioned sources of NG contamination no longer negligible in the future. In other words, the 
same nice properties that make Planck more sensitive to the detection of a primordial NG signal (i.e. high angular 
resolution and sensitivity) make it in fact also much more sensitive to the observation of NG contaminants and require 
a very careful investigation of all the potential sources of bias in the estimate of /nl- 

In this paper we will consider another potential source of uncertainty in the detection of NG, namely the propagation 
of the uncertainties in the cosmological parameters on the measured value of /nl- This effect can be summarised as 
follows: the estimator usually employed to measure /nl assumes a given underlying cosmological model obtained by 
fixing the cosmological parameters at their best-fit values (obtained from the two-point GMB likelihood analysis of 
the experiment under exam). However the cosmological parameter estimates are characterised by uncertainties that 
should be propagated into the /nl estimate in order to accurately quote the final error bars. The uncertainties in 
the cosmological parameters can be safely neglected as long as they are much smaller than the variance of the NG 
estimator. While this works well for WMAP, it is a priori unclear whether it is still a good approximation for Planck. 

The paper is structured as follow: in section [TT] we will describe the /nl estimator commonly employed in the 
analysis and study in detail the effect of cosmological parameters uncertainties on this estimator. After deriving 
the error propagation formulae (|14p and (jl6[) we will apply them to obtain analytical and numerical estimates of 
the expected /nl uncertainty for WMAP and Planck, both in the local and equilateral case. In section Hill we will 
then consider the possibility of applying a more complex analysis in which /nl is estimated by firstly allowing the 
cosmological parameters to vary and then by marginalising over them, rather than by fixing the cosmological model. 
In this case we adopt a Fisher matrix approach to propagate the parameter uncertainties on the final predicted A/nl- 
We will finally discuss our results and draw our conclusion in section ITVl 



II. ERRORS PROPAGATION 



The estimate of /nl from GMB data are usually obtained from measurements of the three point function in harmonic 
space, called the angular bispectrum and defined as: 

KiTr^ ^ • (4) 

Due to rotational invariance of the GMB sky the angular bispectrum can be written as: 

BTXir = orCiTh^.^ia , (5) 

where g^^l':^^' is the Gaunt integral: 



^mim.m3 _ / (2^1 + 1) (2^2 + 1 ) (2^3 + 1) f il h ^3 \ f h h h , 

while bi-^e^i-i is called the reduced bispectrum and contains all the relevant physical information. Analytic formulae 
for both the local and equilateral cases have been computed [H, [lB| . The local reduced bispectrum can be written 
as: 

btti = 2/nT"' / drr'ai, (r)/3,, (r)P,, (r) + (2 perm.) , (7) 
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whereas the equilateral bispectrum is 

bT^iX = S/Z" / drr^^,, (r)/3,, {r)(3,, (r) + (2 perm.) + <5,, {r)S,, {r)5,, (r) + (r)7,, {r)5,, (r) + (5 perm.) . (8) 
The functions ai{r) ^[3i{r) ^^f {r) ^5i{r) appearing in the previous formulae are defined as: 

ai{r) = ^ J dkeAt,{k)jt,{kr) , 

I3i{r) = \ J dkk^P^ik)A,ik)jiikr) , 

7,(r) = ^1 dkePl/\k)A,{k)j,{kr) , 

6,{r) - ^ j dkePl"'{k)A,{k)jdkr) . (9) 

In the previous set of formulae A£(fc) indicates the CMB radiation transfer function and P^{k) is the power spectrum of 
primordial curvature perturbation. It is thus clear that the reduced bispectrum will be dependent on the cosmological 
parameters. 

The NG estimator which is generally employed to analyse CMB data can be written as [13, [ill ■ 

p 1 \ ^ \ ^ /^mimoms ^^i^2^3 /'1n^ 

JNL - 2^ 2^ ^^Ca^s C C C °^i'"i°^2m2Q^3m3 , l^Uj 
^l^2^3 mim2m3 ^3 

where bi-^i^i^ is the analytical form of the primordial reduced bispectrum for the model we are considering (i.e. either 
local or equilateral), whereas is a normalisation factor designed to produce unitary response when /nl — 1: 

^_ (24 + l)(24 + l)(2^3 + l) (li i2 fel')' 

An additional linear term is added to the estimator when dealing with anisotropic noise in the data. We will not 
include such term here as it is not dependent on the cosmological parameters and thus does not affect our estimates. 
The estimated value of /nl is then obtained by correlating the observed bispectrum with the theoretically expected 
one for a given primordial shape (local or equilateral) and given cosmological model, and dividing by a suitable 
normalisation factor. Also the normalisation will be dependent on the bispectrum shape and on the cosmological 
parameters, as both hi^i^i^ and Ci are. The Ci and bi^i^i^ that enter the estimator are calculated by assuming the 
best-fit cosmological model for the experiment under consideration. However, the best-fit cosmological parameters 
are characterised by uncertainties that propagate to /nl- The aim of this work is to analyse this effect in detail and 
assess its significance for WMAP and Planck. As a start, following jT6| let us assume that the cosmological model 
assumed in the calculation of bi^i^i^ and Ci is not the "real" one. Let us call the reduced bispectrum obtained from 
the real cosmological parameters bi-^g^i^. It is then easy to see that the average value of the estimator will be [l6[: 

1 ^ (2£i + 1)(2^2 + 1)(24 + 1) / 4 I2 isV ^W3^(««imia£2m2a£3m3> 



(/nl) - T7 I] 



N ^ 47r VOOO/ Ci,Ct2Ci. 



- y 



{21, + 1)(2^2 + 1)(2£3 + 1) ( h I2 isY fc'^^i^^^a ^^2) 



ei<e2<e3 ^ / «i «2 t3 

throughout the paper we will use a superscript " on /nl whenever we want to indicate a statistical estimate If we now 
want to estimate the bias 5/nl due to the mismatch between the assumed cosmological model and the real one, we 
can write: 

<5/nl= (/nl> -/nl - ^ 1^ ^ [000 C,,Q,Q3 ' 



/nl ^ (2^1 + 1)(2£2 + 1)(24 + 1) fh £2 esV te'te' 



/nl (2^1 + f)(2£2 + 1)(24 + 1) / ^1 h 4 V fe3''^^^i^^^3 . 

tl<l2<l3 ^ / i2 t3 
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In this last formula, Sb^-^i^es is the difference between the theoretical bispectra computed for the "real" and "assumed" 
cosmological model; all the quantities with a superscript ^ are computed in the "real" cosmological model. If we then 
have a cosmological parameter characterised by an uncertainty 6p we simply propagate this uncertainty on the /nl 
estimate as: 



Sp 



6p : 



/nl 
N 



E 

i<i2<i 



(2£i + l)(2£2 + l)(24 + l) 
47r 
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^2 





Sp 



Sp 



(14) 



The term in square brackets expresses the derivative of /nl with respect to the parameter p as a function of the 
derivative of the bispectrum with respect to p. As an estimate of the uncertainty on the parameter p we can use 
the standard deviation and thus substitute 6p, indicating a general small variation of the parameter in the previous 
formula, with tTp, expressing its la uncertainty. The formula above then simply reads: 



a/l 



NL 



dp 



(15) 



this is the standard formula of error propagation. Note that here and in the following equations, we sometimes indicate 
uncertainties in a parameter using the letter S (e.g. Sp), sometimes we use a, e.g. cr/pj^- With S we just generically 
indicate a small variation in the parameter (possibly dependent on the variation in one or more other parameters) , 
whereas with a we specifically refer to the standard deviation. 

In general we have a cosmological model defined by a set of parameters {pi}. All these parameters are allowed to 
vary and can be correlated. The standard error propagation formula in this case becomes: 



\ 



E 



a/] 



NL 



dpi 



dfi 



NL 



Cov(p,;,Pj) 



(16) 



Pj=P3 



where the average values pi of the parameters and their covariance matrix Cov(pj,pj) are determined from a standard 
CMB-likelihood analysis, in which the anisotropies are assumed Gaussian. In our analysis we considered a model 
characterised by the six parameters As,ns,T,uji, = fli,h'^,uj„i — ^mh^ ,^a}, respectively defining the amplitude of 
curvature perturbation at fco =0.002/Mpc, the scalar spectral index, the optical depth to reionization, the physical 
density of baryons, the physical density of matter and the dark energy density. We fixed these parameters at their 
maximum likelihood values from the WMAP 5-year analysis and as an estimate of the covariance matrix we computed 
the Fisher matrix of the parameters. In the computation of the Fisher matrix we consider two cases: in the first case 
we use the 41, 61 and 94 Ghz frequency channels of WMAP (Q+V+W bands) [2lj and in the second we consider the 
combination of the 143 and 217 Ghz frequency channels of Planck [1J|]. For the numerical details of the computation 
(e.g. choice of the step size for the derivatives w.r.t. cosmological parameters) we closely followed the methodology 
presented in We then computed the CMB local and equilateral bispectra numerically using formulae ([7]) , ([5]) , ([5]) 
and we took two sided numerical derivatives to evaluate dB/dpi. For the steps of the derivatives we again followed 
the prescriptions of [l^. Our calculation shows that the /nl error bars get relative corrections ^/^/^'^V.fm^'*' — 16.5% 



and Sf^^ 



cquil / /•cquil 

/ /nl 



14.5% for WMAP in the local and equilateral case, while for Planck we have ^/nl //nl 



,1 / f 1. 

/NL /-/nl 

f local / -flocal 



5% 



and ^/nl"V/nl"' — 4.5%. Before discussing the significance of this correction let us try to understand this result in 
a more intuitive way by employing some analytical approximations. From figures [1] and [5] we see that most of the 
contribution to J/nl from error propagation comes from only 3 of the 6 considered parameters: As^Us and r (this is 
in agreement with the results of [l|]). We will then restrict the following simplified analysis to these three parameters. 

The first step of this analysis is to get analytical expressions for the derivatives of the bispectrum with respect 
to As,ns,T. From formulae (O it can be easily seen that dB/dAs = 2B/As- The fractional variation in /nl 
corresponding to a variation SAs is then: 



Shh . ■ „ ; SAs 
-SAs = 2/nl- 



SAs 



As 



(17) 



The parameter t defines the optical depth to reionization and the effect of changing it can be described by introducing 
a multiplicative factor e""^ in front of the radiation transfer function at high £'s. Note that the radiation transfer 
functions appear in the definition of the bispectrum through the functions a, (3, 7, S defined in formula ©. The 
bispectrum modes that give the largest contributions in the local case to the final signal-to-noise ratio are the so 
called squeezed configuration, i.e. configuration where one of the three £'s is much smaller than the other two, and 
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FIG. 1: Contribution of the different cosmological parameter uncertainties to the final error in the estimate of /nl^' (lower 
panel) and /J^l"' (upper panel) as a function of ^max- The quantity (S/nl/^pO""?; plotted for each of the six parameter in 
the model. In this figure we considered an experiment with the characteristics of WMAP. 



equilateral triangles in the other case (see e.g. [16]). For this reason, in the local case one of the three £'s will be super- 
horizon and the corresponding transfer function will not show a multiplicative factor e^'^ in front. In the equilateral 
case all modes are sub-horizon in the important configurations. We can then write b^gl^gf^^ = exp(— 2t)&^°'^^'^^ and 
6°^;^^'^^ = exp(-3T)6^;i;;i^. Substituting into equation ^ we find Sfj^l""^ ~ -2fj^l''^5T and (5/^1"'' -Sf^'^j^'^ST. 
Note though that the parameter t obtained from the Cg likelihood analysis is degenerate with the amplitude of the 
spectrum of primordial curvature perturbations. In order to include this degeneration in our simplified description, 
for a given variation in r we will also introduce a variation in the power spectrum amplitude that leaves the final 
Ci unchanged. This is obtained by multiplying the amplitude by a factor exp(a(5r), where a = 2,3 in the loci and 
equilateral case respectively. A small shift St in the ionisation optical depth then implies a shift in the amplitude 
equal to 6 As — 2AsSt. The total bispectrum variation SB is then given by (we omit the subscript ^i£2^3 for simplicity 
of notation): 

5B ^ dB ^ dB 5 As , 
or OT oAs OT 

= -aBe-^^^Sr + ABe''^^ST 

~ {A-a)B6T, (18) 
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FIG. 2: Contribution of the different cosmological parameter uncertainties to the final error in the estimate of /nl^' (lower 
panel) and /^l"' for an experiment with the characteristics of Planck. 



where in the last line we neglected second order terms in St. The total variation in /nl for a given 5t is then: 



St 



-,5t= (4-a)/NL<5T (19) 



The remaining parameter to take into account is the scalar spectral index n. Our next step is then the evaluation of 
SB/Sn. First of all we note that when changing n we have to change the power spectrum normalisation accordingly 
because the normalisation is defined at a given pivot scale. To compute SB arising from a small change in the spectral 
index we then have to evaluate again: 

SB^ dB ^ dBSA^ 

where the partial derivative with respect to n is taken by assuming A fixed. The authors of [13] use WMAP 3-years 
data to find that the normalisation is well fit by the following expression: 

j^WMAP ^ ^ exp(~1.24+l.G4r)(l-n) 
^ ^ Vl + 0.53r ' ^ ' 
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where r is the tensor-to-scalar ratio. We wiU use this ansatz, with the additional assumption r = 0. In this way we 
obtain, for a given variation 6n of the scalar spectral index: SA = 1.24:A6n and, correspondingly, 6B ~ 2.5dn. Finally, 
to approximately evaluate dB/dn we work in the pure SW regime. Estimates of the signal-to- noise ratio have in this 
case been obtained by Komatsu and Spcrgel for the equilateral configurations and by Babich and Zaldarriaga [Tsl 
for the squeezed ones for n = 1. Extending their results, we obtain that in both cases: 



NL 



dn 



2 



l0g(^n 



1 



(l-n) 



(22) 



this allows to write the variation in /nl for a given 6n 



/nl 



2-5 - 7777^ + ^ log (^max) 

2(1 — n) 2 



Sn . 



(23) 



Having an expression for the derivatives of the bispectrum with respect to each of the three parameters As,n and r 
we can now propagate the error using Eq. (|16|). that for this particular case reads 



a; 

JNL 



NL 



^1 

-45 



NL 



6t 



NL 



5n 



(24) 



Note that, as mentioned above, the correlation of As with r and n is accounted for by the ansatz we have made for 
As: 



As = As exp [-ar - 1.24(1 - n)] 



(25) 



whereas the correlation between r and n has been neglected in the previous formula. Using the expressions just 
derived above for (5/nl/^^S: ^Inl/^t and Sfj^i^/Sn we finally get: 



JNL 



/l 



NL 



2^) +[{^~a)arf 



2.5 - 
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2(1 



^fog( 







(26) 



The present WMAP 5-years analysis yields a fractional uncertainty on the amplitude of the curvature power spectrum 
of order 3%, while cr,- ~ 0.016 and ct„ = 0.015. Substituting these numbers in the last formula yields a total fractional 
correction of order 14% on /nl, in very good agreement with the numerical results. If we now consider an experiment 
with the characteristics of Planck, our Fisher matrix analysis give a fractional uncertainty on ^5 of order 1.5%, 
cr„ = 0.004, ar — 0.005. This produces a final fractional correction of order ~ 5% on /nl, again in very good 
agreement with the numerical estimate. 

To understand whether these corrections are negligible or not we have now to compare it with the la uncertainty 
^/nl of the estimator, obtained with fixed cosmological parameters. The WMAP analysis finds A/^^'^' ~ 30, 
^/nl"' — 100, and central values /nl''' — 60 and /nl"'' — '''0. Using the fractional uncertainties above we get a 
contribution to the final error bar from cosmological parameters uncertainties that amounts to A/jlf^^' ~ 10 and 
^/nl"' — 9- The effect of propagating cosmological parameters uncertainties can then be neglected for the equilateral 
shape, where the error bars are larger, whereas it is more important for the local shape (about 30% of the presently 



quoted error bar). If we now consider Planck, Fisher matrix based forecasts predict A/ 



local 
NL 



5 and A/ 



cquil 
NL 



60. 

If we assume the /nl central values found by the WMAP analysis we obtain a small effect for the equilateral case, 
whereas 5f^£'^^ ~ 3, i.e. of the same order of magnitude as Af^f^^^. This analysis then suggests that the effect of 
propagating uncertainties in the cosmological parameters on the final /nl**' error bar should be taken into account 
if large central values of /n^ ' are found with Planck. Note that a value of /nl'^' of order 60 would mean a many a 
detection with Planck. Correcting the error bar in order to account for error propagation effects would not change 
this result but it would on the other hand modify the level of significance of such a detection. 

Before concluding this section, we would like to stress again that the estimator of /nl currently employed in the 
analyses fixes the cosmological parameters at their best-fit values. A way to reduce the impact of the uncertainties 
on the parameters would be to perform a joint likelihood analysis in which the cosmological parameters are allowed 
to vary and then marginalise over their uncertainties. Obtaining a forecast of the final /nl error if this approach is 
taken is the purpose of the next section. 
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III. FISHER MATRIX 



As we were mentioning in the previous section, the optimal approach to the /nl measurement would be to treat 
the cosmological parameters as nuisance parameters and to marginalise over their distributions in order to get the 
final /nl estimate. The error on /nl can in this case be estimated by a Fisher matrix analysis. If we consider a set 
P = {Pi} of cosmological parameters we can express the Fisher matrix as p^. [23|: 

~ ^ dp, dp, ' ^^'> 

where is the bispectrum variance. In the limit of small non-Gaussianity we can take — C i^C g^C g^^i^i^g^ , where A 
takes the values 1,2,6 when two £'s are different, two of them are the same and all are the same respectively. Following 
the results of the previous section, we know that the relevant set of parameters to consider is p = {/nl, 't-, t}. 
Two account for cosmological parameter uncertainties we add a Gaussian prior on the ith parameter with variance 
erf, where ai is the standard deviation obtained from the two-point function likelihood analysis. This approach is 
feasible as long as we deal with weak non-Gaussianity and the two and three point function can then be treated as 
uncorrelated. A Gaussian prior on the ith parameter with variance of is imposed by simply adding a ^ term to 

the a entry of the Fisher matrix (see e.g. [Ill)- Once the Fisher matrix has been computed, the error on the i-th 
parameter after marginalising over the others can be estimated in the standard way as: 



^7' ■ (28) 

Before moving to the numerical evaluation of formula (|27p for the full set of parameters let us start with a simplified 
case in which only /nl and r are considered in the analysis and let us for simplicity restrict ourselves to the local 
case. In this case, having made the approximation (explained in the previous section) dB /dr = exp(— 2t)_B, a simple 
analytical calculation gives the following Fisher matrix: 



F = ^ 2 „2 , (29) 

9 f.T. -"^l'^2<!3 Af2 -"^1'?2<'3 

-^./nl 2^ ^-i^ 4/nl „i 

where ct^ is the bispectrum variance defined above. This matrix is singular for /nl 0, meaning that /nl and r are 
degenerate parameters. Adding a Gaussian prior on r with variance cr^ breaks the degeneracy. 




^ 2 „T . (30) 

V-2/nlE^ 4/5LE^ + ^y 

Inverting the Fisher matrix and taking the square roots yields the final error on /nl: 

= , TTl^ \ i + 4f7^/^L } ] ■ (31) 



If we call A/nl the estimated Fisher matrix error when we do not marginalise over t (i.e. the error usually quoted in 
the literature) then we see from the previous formula that: 



y'(A/NL)2 + 4a2/2^; (32) 



note the difference with respect to the previous approach in which cosmological parameters were fixed. In that case 
the cosmological parameter errors biased the estimator and the uncertainties propagated linearly (see also p^): 

^/nl = ^/nL + 2/NLO-r ■ (33) 

As we saw in the previous section, the error propagation scheme arising from the standard approach of fixing cos- 
mological parameters produces a relative correction of a few percent for WMAP and Planck. We concluded that 
this correction is small but not always negligible for Planck. On the other hand the marginalisation approach used 



9 




FIG. 3: Correction to A/nl after marginalisation over r is performed in a toy model where only uncertainties over r are 
considered. We plot the correction as a function of A/nl. The correction becomes significant only when A/nl is small enough 
to produce a many g detection for a given /nl- 



here makes the additional uncertainty much smaller than it was previously and always negligible, even for very large 
central value of /nl, unless /nl is large enough to produce a many a detection for a given experiment (see Fig. [S]). 

As long as only the parameter r is considered we can then conclude that both for WMAP (ctt = 0.016, A/^^**' ~ 30) 
and Planck ((Tt — 0.016, A/^^^' ~ 5) the correction to the /nl error bars is totally neghgible: Sf^l''^ < 0.2%/nl, 
assuming a central value /nl — 60. From the formula above we basically see that the effect of marginalising over 
r is to suppress the correction mentioned in the previous section by a further factor A/nl- Moreover we recover 
the correction mentioned in the previous section in the limit A/nl 0. All this makes sense: the correction from 
cosmological parameters uncertainties is significant only if the error bar on /nl before marginalisation is comparable 
to the error bars on the other parameters; moreover a full likelihood estimation optimises the final error bar on /nl 
with respect to an analysis in which the cosmological parameters are held fixed. The same results arise when we 
account not only for r, but we consider the full set {^,n,r, /nl}. In this case we evaluated (5/nl numerically from 
formula (|27p and obtained that the correction on the /nl error bar after marginalisation is always less than 0.5%. 
The conclusion is that if a full likelihood analysis including the two and three point functions is applied in order to 
estimate /nl, then the impact of cosmological parameters uncertainties is totally negligible. 
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IV. CONCLUSIONS 



In this paper we considered the effect of propagating cosmological parameters uncertainties on the estimate of the 
primordial NG parameter /nl- We firstly show that, accounting for the large central value of /nl presently measured 
[l|, Q , the final correction from parameters uncertainties is of order 30% of the quoted f^{f^ error bar for WMAP and 
at about the same level of the predicted /^l''' error bars for Planck. If a large /nl will be observed by Planck, the 
effect of these uncertainties will then be not big enough to change the conclusion that a large level of primordial non- 
Gaussianity is present in the data. However the effect is important enough to change the significance of the detection 
and should be taken into account when quoting the error bars. We finally show that the effect of cosmological 
parameters uncertainties becomes totally negligible if we do not fix the cosmological parameters in the analysis, but 
we treat them as nuisance parameters and marginalise over their distribution in order to obtain the final /nl estimate. 
Even if optimal, this last approach is nevertheless probably still inconvenient. A joint-likelihood evaluation would 
require a large amount of time and the final gain in the error bar would be significant only for large values of /nl, 
but those would produce a significant detection even in the sub-optimal approach. 
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